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Abstract 

Conditions during fetal development influence health and disease in adulthood, especially during critical windows of 
organogenesis. Fetal exposure to the endocrine disrupting chemical, bisphenol A (BPA) affects the development of multiple 
organ systems in rodents and monkeys. However, effects of BPA exposure on cardiac development have not been assessed. 
With evidence that maternal BPA is transplacentally delivered to the developing fetus, it becomes imperative to examine 
the physiological consequences of gestational exposure during primate development. Herein, we evaluate the effects of 
daily, oral BPA exposure of pregnant rhesus monkeys {Macaca mulatto) on the fetal heart transcriptome. Pregnant monkeys 
were given daily oral doses (400 |xg/kg body weight) of BPA during early (50-100±2 days post conception, dpc) or late 
(100±2 dpc - term), gestation. At the end of treatment, fetal heart tissues were collected and chamber specific 
transcriptome expression was assessed using genome-wide microarray. Quantitative real-time PCR was conducted on select 
genes and ventricular tissue glycogen content was quantified. Our results show that BPA exposure alters transcription of 
genes that are recognized for their role in cardiac pathophysiologies. Importantly, myosin heavy chain, cardiac isoform 
alpha (Myh6) was down-regulated in the left ventricle, and 'A Disintegrin and Metalloprotease 12', long isoform {Aclaml2-I) 
was up-regulated in both ventricles, and the right atrium of the heart in BPA exposed fetuses. BPA induced alteration of 
these genes supports the hypothesis that exposure to BPA during fetal development may impact cardiovascular fitness. Our 
results intensify concerns about the role of BPA in the genesis of human metabolic and cardiovascular diseases. 
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Introduction 

Bisphenol A (BPA) is used in tiie manufacture of polycarbonate 
plastics and epoxy resins that line food and beverage containers 
and a wide variety of consumer products such as water bottles and 
sales receipts; this results in repetitive, long-term exposure in most 
humans. BPA exposure in pregnant women can cross transpla- 
centally, potentially impacting the developing fetus [1,2]. While 
the adult conjugates BPA into BPA-glucuronide (BPA-GA) in the 
hver for excretion in the bile, this metabolic pathway is quite low 
in the fetus; of greater concern, BPA-GA is deconjugated in the 
fetus increasing exposure to the more biologically potent BPA [2] . 
Furthermore, BPA can interact with several biologically important 
hormone receptors; estrogen [3,4], androgen [5], thyroid [6,7] and 
intercellular communication channels such as CX-43 based gap 
junctions [8], though it remains unclear whether the associated 
signaling cascades are activated or inhibited. There is growing 
evidence that the receptor-mediated responses to BPA are 
dependent upon dose, duration and timing of BPA exposure as 
well as the cell/tissue studied [9-12]. 



The Barker hypothesis, first presented to explain the increased 
risk of cardiovascular disease in children of undernourished 
mothers of the Dutch famine, suggests that fetal developmental 
"programming" predisposes some individuals to the development 
of chronic disease later in life. The crux of this theory is that the 
phenotype of an individual is dependent upon both genotype and 
environmental influences that alter prenatal growth rate, temporal 
growth patterns and birth weight [13]. Most relevant to this work 
is that Barker extended this theory by directly linking the incidence 
of cardiovascular disease and insulin resistance to the placental 
environment [14]. While the BPA specific effects on adults is 
widely acknowledged in animal studies [15-17], fetal BPA 
exposure also affects multiple organ systems, with behavioral, 
reproductive and metabolic consequences in the adult [18]. 
Recent evidence also show that exposure to BPA during 
pregnancy in mice not only alters glucose homeostasis but also 
serves as a risk factor for subsequent development of diabetes and 
heart disease in the offspring [19]. Intriguingly, epidemiological 
studies also suggest a correlation between higher urinary BPA 
concentrations in adult men and increased incidence of cardio- 
vascular disease and diabetes [20-22] 
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Cardiovascular disease and diabetes are chronic Ulnesses tliat 
are rising in incidence and are among the top 5 diseases globally in 
terms of burden on the health care system. Cardiac tissue is known 
to express estrogen [23], androgen [24] and thyroid receptors 
[25,26] and these hormones play critical roles in the normal 
developmental process of cardiac tissue; as such we investigate the 
potential for disruption of these developmental mediators by 
exposure to BPA. The objective of this study was to evaluate the 
eflFects of daily, oral maternal BPA exposure in primates on the 
gene expression profiles in each of the left ventricle (LV), right 
vcntrick- (RY), left atrium (LA) and right atrium (RA) of 
developing fetal heart tissue. This study is the first report of the 
impact of maternal BPA exposure on fetal/neonatal cardiovascu- 
lar tissue in a primate. 

Materials and Methods 

Ethics statement 

All animal procedures were approved in advance by the 
Institutional Animal Use and Care Advisory Committee at the 
University of Cahfomia, Davis, and all procedures complied with 
the requirements of the Animal Welfare Act. 

Animal experimentation and tissue collection 

The tissues evaluated in this study were part of a study to 
determine the effects of BPA on ovarian development and detailed 
information about animal care and BPA exposure have been 
reported [27]. Therefore, all analysis reportc'd in this study were 
from female fetal heart tissue. All animal procedures were 
approved in advance by the Institutional Animal Use and Care 
Advisory Committee at the University of California, Davis and all 
procedures complied with the requirements of the Animal Welfare 
Act. Adult female rhesus monkeys (Macaca mulatto) were housed at 
California National Primate Research Center and maintained 
according to the institutional guidelines (http://www.cnprc. 
ucdavis.edu) as previously described [27-29]. Animals were 
housed individually in stainless steel cages and maintained under 
a 0600-1800 hours light cycle at 25-27°C. All animals were 
allowed ad libitum acc(;ss to Purina Monkey Chow and water 
supplied through PVC pipes and a "lixit" device. Further, 
environmental enrichment was evoked by offering supplemental 
seasonal produce, seeds and cereals. Animals with a history of 
prior pregnancy were bred using established time-mated proce- 
dures and pregnancy was detected by ultrasound [30,31]. 
Pregnant monkeys were sonographically screened early in 
gestation to confirm normal development and to identify female 
fetuses (by 40 days gestation) using estabhshed protocols with 
proven reliability [30,32]. Fetuses were delivered vaginally or 
through a Caesarian section. Fetal tissues were collected at 
necropsy and euthanasia conformed to the recommendations of 
the American Veterinary Medical Association (AVMA) Guidelines 
on Euthanasia and Primate Center SOPs (overdose of pentobar- 
bital). 

The BPA dose selected for this study was determined to lead to 
peak blood levels in non-pregnant monkeys [28] that approximate 
those levels that are commonly measured in humans [33]. 
Pregnant females were administered a small piece of fruit once 
daily that contained either deuterated BPA (400 Hg/kg body 
weight), dissolved in 150 ^ll ethanol (vehicle) or 150 |J,1 ethanol 
only (control), during either early gestation, EG (days; 50-100±2; 
n = 5) or late gestation, LG (days; 100±2 - term; control, n = 5 and 
BPA, n = 4) [28] . Quantification of maternal serum BPA indicate 
peak levels of conjugated BPA of 14 1.08 ±26.64 ng dBPA/ml and 
unconjugated BPA of 1 .99±0. 1 2 ng dBPA/ml in the EG timeline. 



In the LG time line, the maternal serum conjugated and 
unconjugated BPA levels were 39.09±15.71 ng/ml and 
0.68±0.312 ng/ml, respectively [29]. As explained earlier, fetal 
heart tissues of left ventricle (L\'), right ventricle (R\'), left atrium 
(LA) and right atrium (RA) were obtained from fetuses removed by 
Caesarian section (EG; control, n = 5 & BPA, n = 4) or ax necropsy 
following vaginal delivery (LG; control, n—5&. BPA, n — 4). Other 
fetal tissues were harvested from these fetuses and results of their 
analyses have been reported previously [27,29]. 

RNA isolation 

Total RNA was extracted from rhesus monkey fetal LV, RV, LA 
and RA tissue using TRIzol reagent (Imdtrogen, Carslbad, CA, 
USA) according to the manufacturer's protocol. RNA was 
quantified using Nanodrop® ND- 1 000 UV-Vis Spectrophotometer 
(Nanodrop Technologies, Wilmington, DE, USA) and integrity of 
RNA was assessed using 1 .5% denaturing formaldehyde agarose gel 
electrophoresis. Subsequently, 15 ^ig of DNase treated (Ambion, 
Foster City, CA, USA) total RNA was utilized for cDNA synthesis. 

Custom NimbleGen Rhesus Microarray design 

A custom unbiased NimbleGen rhesus macaque 1 2-plex whole 
transcript gene expression microarray (Roche NimbleGen, 
Madison, WI, USA) was designed from the European Bioinfor- 
matics Institute (FBI) Ensembl database [Macaca mulatta, Mmul_l, 
Jan 2006). We identified gene transcripts (21,905) and small RNAs 
(5,253) from the entire genome (total 28,920 potential transcripts). 
Of the 21,905 transcripts targeted, we were able to build probes 
for 21,840 of which 21,808 contained 6 probes each and the rest 
(32) had less than 6 probes. Of 5,253 sRNA targeted, we were able 
to build probes for 4,309 of which 1,974 had a single probe and 
2,335 contained 2 probes for a total of 137,577 unique probes (60 
'mer in length) on the array. 

Microarray sample preparation, hybridization and 
scanning 

Briefly, 10 ng of DNase treated total RNA from LV or RV was 
primed with 100 picomoles of oligo (dTjs) and 200 ng of random 
primers and reverse transcribed using 1.5 |iL of superscript III 
(Invitrogen). Second strand cDNA was synthesized and RNase 
treated to remove residual RNA. Due to the limitations in the 
quantity of tissue available from the LA and RA, cDNA 
preparations from both these tissue were prepared by employing 
a PGR based amplification approach available through Trans- 
Plex® Complete Whole Transcriptome Amplification Kit (WTA2) 
(Sigma, St. Louis, MO, USA). In this case, double stranded cDNA 
(dscDNA) was synthesized from a 40 ng of DNase treated total 
RNA from each sample. In both cases, resultant dscDNA was 
purified by employing phenol:chloroform:isoamyl alcohol, ethanol 
precipitation and quantified. Finally, cDNA was labeled using Cy5 
nanomers (NimbleGen). One 1 2-plex array that contained 12 sub- 
arrays was used per tissue (LV, RV, LA & RA) to hybridize three 
samples (control, n=3 & BPA, n = 3) from each of the BPA or 
control treatment groups, in both the EG and LG time lines (2 
(LG, EG) X 2 (control, BPA) x 3 (replicates) = 12 samples). Labeled 
cDNA (2 ng) from each sample was hybridized overnight at 42°C 
using the NimbleGen Maui Hybridization .system 4 (NimbleGen). 
Arrays were washed using a Maui Wash system (BioMicro 
Systems, Salt Lake city, UT, USA). Genepix pro 6.1 software 
(Molecular devices, Downingtown, PA, USA) operating a 
Genepix4000B scanner (Molecular devices) was used to scan all 
hybridized arrays. Scanning was performed at a resolution of 
5 nm, and photomultipHer tube (PMT) gains were set to obtain 



PLOS ONE I www.plosone.org 



2 



February 2014 | Volume 9 | Issue 2 | e89096 



Maternal BPA Exposure Alters Fetal Cardiac Genes 



maximum signal intensities such that approximately 5% of the 
probes reached signal saturation (2"^ raw signal intensity). 

Quantitative real-time Polymerase chain reaction (qRT- 
PCR) 

DNase treated total RNA (2 |ag) from each tissue sample of both 
EG (control, n=5 &. BPA, n = 5) and LG (control, n = 5 and BPA, 
n = 4) timeline was reverse transcribed using high-capacity cDNA 
reverse transcription kit (Applied Biosystems, Foster City, CA, 
USA). Gene expression was measured by employing custom 
designed gene specific Taqman MGB® primer/probes on an ABI 
7500 fast real-time PGR system (Applied Biosystems). Primers 
(Invitrogen) and probes (Applied Biosystems) (Table 1) were 
designed to span more than one exon of a gene, using Primer 
Express 3.0 (Applied Biosystems). Real-time PGR assays were 
conducted as reported previously [34]. AU genes were assayed in 
duplicate. Endogenous reference, 18 S rRNA was measured using 
0.5 ng cDNA, and all other genes were assayed using 50 ng cDNA 
in each reaction. 

Glycogen analysis 

Glycogen content in the ventricle (LV, RV) of the LG group 
fetuses (control, n = 5 and BPA, n = 4) was assessed using a 
modified colorimetric method [35] . Pre-weighed ground tissue (5- 
20 mg) in 0.3 ml of 30% KOH was boiled for 30 minutes with 
mfadng. Samples were cooled to room temperature, 0.1 ml of 1 M 
Na2S04, and 0.8 ml of 100% ethanol was added and samples 



Table 1. Taqman MGB primer & probe system. 









Primer & Taqman Probe 


Gene 


Gene accession number 


sequences 


1 8 s rRNA 


AF243428 


FP: CCACGCGAGATTGAGCAAT 


RP: GCAGCCCCGGACATCTAA 


TP: ACAGGTCTGTGATGCC 


Myh6 


ENSMMUT00000020702 


FP: AGCGGCGCATCAAGGA 






RP: 






CGCAGCAGGTTCTTTTTGTCT 


TP: CTCACCTACCAGACAGAG 


Adam 12-1 


ENSMMUT00000013487 


FP: CGGCCACCTCGGAAAAG 






RP: 






CCTTCGGTGGGTAGGAATCTG 


TP: CCTGATGAGGAAGCC 


{Thrb) 


ENSMMUT00000000096 


FP: TGAATCAATACAGCCAAC 






CTGAA 






RP: GGCGACTGCACTTGAGAAA 






AA 


TP: ATTTCACTGAAGAAAAGC 


Gygl 


ENSMMUT00000007776 


FP: CTACTCCTCAGGTCTCAGA 






CTCCAT 




ENSMMUT00000007779 


RP: CTGTCCAAGACATCTACCAT 






GATGA 




ENSMMUT00000007773 


TP: AAGTTTTAGAGACAGTCTTC 


Pdk4 


ENSMMUT00000030108 


FP: TGCACCAACGCCTGTGAT 






RP: GCAAGCCGTAACCAAAAC 






CA 


TP:TTCCCGGAATGCTCCT 


Sequences of primers; forward (FP), reverse (RP), Taqman probe (TP) specific to 



gene transcripts analyzed. 
doi:l 0.1 371/journal.pone.0089096.t001 



were boiled for 5 additional minutes. Samples were centrifuged at 
15000 X g for 10 minutes at 4°G, supernatant was removed and 
pellets were washed twice in 0.2 ml of MiUipore water and 0.4 ml 
of 100% ethanol. Pellets were resuspended in 0.2 ml of 
amyloglucosidase solution (0.3 mg/ml amyloglucosidase in 
0.2 M NaOAc pH 4.8), incubated at 40°C for 3 hours. Subse- 
quently, the reaction mix was diluted with 0.2 M NaoAc; pH 4.8. 
Quantification of glucose, derived using hexokinase and glucose 6 
- phosphate dehydrogenase from tissue glycogen, was measured at 
340 nm using Evolution 300 UV -Vis spectrophotometer (Thermo 
scientific, Waltham, MA, USA). 

Microarray data analysis 

NimbleScan software (NimbleGen. Madison, WI) was used to 
align a chip-specific grid to control features and extract raw 
intensity data for each probe and each array. Raw intensity data 
were then read into the R statistical computing environment and 
checked for quality (R Development Gore Team [2010]: A 
language and environment for statistical computing R Foundation 
for Statistical Computing, Vienna, Austria ISBN 3-900051-07-0 
URL http://www.R-project.org). Further, chip intensity distribu- 
tions, boxplots and hierarchical clusters were compared and 
evaluated for any unusual global patterns. Within each tissue, each 
array was background corrected and normalized using the 
quantile normalization procedure and each probeset was summa- 
rized using the median pohsh procedure as described by the robust 
multichip average (RMA) procedure. Probesets with low levels of 
expression variation across all samples (IQR<0.5) were removed 
from further analysis, reducing the overall number of statistical 
tests to be performed. Differential expression was assessed using a 
linear model with an empirical Bayesian adjustment to the 
variances and comparisons of interest were extracted using 
contrasts (n = 3). Genes that differ by an arbitrary differential 
expression filter criteria (DEFC) of logj fold change (LEG) = ± 1 at 
an unadjusted ^^0.01, were considered significantly altered in 
expression. Both the raw and RMA normalized microarray data 
were deposited in the Gene Expression Omnibus (GEO) database 
(accession number. GSE53393). 

Quantitative real-time PCR and glycogen data analysis 

Taqman based gene expression analysis using qRT-PGR was 
performed as described previously [34,36]. Briefly, delta Gt (AGt) 
values for each gene and aU samples were obtained through 
subtracting the matched 18 s rRNA values from the averaged 
threshold cycle (Gt) values of candidate genes. Delta Gt values 
(AGt) served as the response variables (control, n = 5 and BPA, 
n = 4) for statistical comparisons. DifiFerential average response of 
each gene for BPA exposed over control (AAGt) is reported as 
mean log2 fold change (LEG) ± SEM, and represented graphically 
using Sigma plot 8.0 software (SPSS Inc., Chicago, IE, USA), 
Tissue glycogen content was expressed as |j.g/ mg wet tissue weight. 
Data were analyzed using SAS software (Statistical Analysis 
Software, Gary, NG, USA). A two sample t-test was performed to 
identify mean significant differences in the qRT-PGR measured 
gene expression and the glycogen content between the BPA 
exposed fetuses and the matching control group. A two-tailed 
/)^0.05 was considered statistically significant. 

Results 

Gene expression changes observed using whole genome 
DNA microarray 

Microarray analysis of cardiac tissues (LV, RV, LA & RA) 
indicated a differential mRNA expression profile in the EG and 
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the LG maternal BPA exposed fetuses compared to that of 
matching controls. Genes that changed by &2 fold (LFC = ±1) at 
p^O.Ol (unadjusted /^-m/mi?) were considered significantly altered 
by maternal BPA exposure. This resulted in lists of up- or down- 
regulated genes in cardiac tissues analyzed after EG (Fig. lA; 
Table S1-S4) or LG (Fig. IB; Table S5-S8) maternal BPA 
treatment. Differentially expressed genes included both protein 
coding genes and a variety of non-protein coding RNAs; micro 
RNAs (miRNAs), small nucleolar RNAs and other RNAs. Manual 
screening of the difiFerentially expressed gene lists identified two 
important protein coding genes that changed significantly; Myh6 
and Adaml2, and both these genes were recognized for their role in 
cardiac pathophysiologies. Our results show that Myh6 is 
significantiy down-regulated in the fetal LV (LFC = —5.385, 
—42 fold), but not in the RV or the atria, after LG maternal 
BPA exposure, or in any tissue type after EG maternal BPA 
exposure, compared to the matched controls (Table 2-3). On the 
other hand, Adaml2 is significantly up-regulated in the LV 
(LFC = 1.678, «3.2 fold) and the RV (LFC = 1.213, «2.3 fold) 
of the LG, but not the atria of the LG, or in any tissue of the EG, 
maternal BPA exposed fetuses compared to matched controls 
(Table 2-3). 

Thyroid hormone status and glycolytic flux are important 
determinants of Mjh6 expression (Gupta 2007). Accordingly, 
thyroid hormone receptors play an important role in the Mjh6 
expression, and thyroid hormone receptor beta 
(Q8SPN6_MACMU - Thrb) is up-regulated in the LV of the 
LG, BPA treated fetuses vs. matched controls (Table 3). However, 
this did not meet our differential expression filter criteria as the 
unadjusted p = 0. 055. Also, preliminary studies in our lab indicate 
that maternal BPA exposure may interfere with the glycogen 
metabolism in mice fetal heart tissue as measured through altered 
pyruvate dehydrogenase kinase isozyme 4 {Pdk4) and glycogenin 
{Gygl) mRNA expression (data not shown). Microarray analyses 
revealed no significant differences in the mRNA expression levels 
of either Pdk4 or Gygl between the maternal BPA exposed v.s. 
matched control fetuses (Table 2-3). 

We also observed induction of two miRNAs; miR-205 
(LFC = 2.583, =6 fold) and miR-224 (LFC = 4.093, «17 fold) in 
LV of the LG, maternal BPA exposed fetuses compared to that of 
the matching controls (Table 2-3). Collectively, microarray results 
indicate that maternal BPA exposure alters the transcriptional 



Table 2. Differentially expressed genes observed using 
microarray on rhesus fetal heart exposed maternally to BPA or 
vehicle during early gestation, EG (days; 50-1 00±2). 





CMCCiuiRL in 


Gene 


EG - loc 


2 fold change 








LV 


RV 


LA 


RA 


ENSMMUG00000014754 


Myh6 


-0.204 


2.619 


ND 


ND 


ENSMMUG00000009661 


AdamU 


0.09 


-0.115 


-0.225 


0.192 


ENSMMUG00000000067 


Thrb 


0.845 


-0.441 


ND 


0.536 


ENSMMUG00000005521 


Gygl 


ND 


-0.177 


ND 


ND 


ENSMMUG00000021389 


Pdk4 


0.498 


0.188 


0.016 


-0.263 


ENSMMUG00000026923 


mmi-mir-205 


-0.183 


ND 


-0.385 


0.403 


ENSMMUG00000026877 


mm!-mir-224 


0.581 


ND 


-0.459 


-0.069 



ND - below differential expression filter criteria. 
doi:1 0.1 371 /journal.pone.0089096.t002 



Table 3. Differentially expressed genes observed using 
microarray on rhesus fetal heart exposed maternally to BPA or 
vehicle during early gestation, LG (days; 100±2-term). 



ENSEMBLJD 


Gene 
symbol 


LG - log2 fold change 








LV 


RV 


LA 


RA 


ENSMMUG00000014754 


Myh6 


-5.385 


-0.187 


ND 


ND 


ENSMMUG00000009661 


AdamU 


1.678 


1.213 


0.091 


0.555 


ENSMMUG00000000067 


Thrb 


1.185 


-0.865 


ND 


-0.096 


ENSMMUG00000005521 


Gygl 


ND 


-0.143 


ND 


ND 


ENSMMUG00000021389 


Pdk4 


0.346 


0.171 


0.22 


0.322 


ENSMMUG00000026923 


mml-mir-205 


2.583 


ND 


-0.909 


0.109 


ENSMMUG00000026877 


mml-mir-224 


4.093 


ND 


0.275 


0.161 



ND - below differential expression filter criteria. 
doi:1 0.1 371/journal.pone.0089096.t003 
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Figure 1. Total number of significantly altered genes observed using microarray global transcriptome expression analysis. Whole 
transcriptome analysis was performed on the left ventricle (LV), right ventricle (RV), left atrium (LA) and the right atrium (RA) of rhesus monkey 
(Macaca mulatta) fetuses that were exposed maternally to a 400 ng/kg body weight, Bisphenol A (BPA), relative to matched control fetuses, either 
during, (A) early gestation, EG (days; 50-100±2) or, (B) late gestation, LG (days 100±2-term). Genes that changed by greater than ±1 log2 fold 
change (1 LFC = 2 fold change) at an unadjusted p<0.01 were considered significant and differentially expressed. Bars represent the number of 
upregulated (positive y-axis) or downregulated (negative y-axis) genes in each of the tissues. 
doi:1 0.1 371/journal.pone.0089096.g001 



PLOS ONE I www.plosone.org 



4 



February 2014 | Volume 9 | Issue 2 | e89096 



Maternal BPA Exposure Alters Fetal Cardiac Genes 



profile in rhesus fetal heart pertaining to cardiac remodeling events 
that have been linked to cardiovascular pathophysiologies. 

Quantitative real-time PCR analysis 

Using qRT-PCR, we quantified the transcript abundance of 
selected genes in all the heart tissue types from the LG maternal 
BPA exposed fetuses (n = 4) and the matching controls [n = 5) 
(Table 1). Myosin isoform Mjh6 specific mRNA expression was 
measured. Two isoforms oiAdam.12, the membrane anchored long 
isoform {Adaml2-l) and the secreted short isoform {Adaml2-s), arise 
through alternative splicing ot Adam 12 transcript. Using a specific 
primer probe set for Adaml2-l (Table 1), we measured and report 
the transcript abundance oi Adaml2-l isoform. Similar results were 
obtained when we used another primer-probe set that simulta- 
neously detects both isoforms (data not shown) indicating the 
predominant expression of the long isoform in these tissues. 
Although the Thrb, Gygl and Pdk4 genes did not pass our filtering 
or significance criteria in the microarray study (Table 2-3), we 
elected to quantify mRNA expression for these genes using the 
more sensitive qRT-PCR. Expression of Af))/! 6' mRNA was down- 
regulated by «38 fold (LFC = - 5.243±2.544, = 0.077) in the 
LV of maternal BPA exposed fetuses compared to matching 
controls (Fig. 2). While the difference is marginally significant, we 
believe this observation is important to report since the variability 
is largely caused by one fetus in each group. Also, Adaml2-l 
expression increased similarly in both ventricles, reaching signif- 
icance in the RV (LFC= 1. 528±0. 6039, =2.9 fold, /) = 0. 031) 
and trending in the LV (LFC = 1.57± 1. 1087, «3 fold, = 0.154), 
of the BPA exposure group vs. matched control (Fig. 3). Increased 
Adaml2-l expression was also observed in the RA of the BPA 
treatment group compared to controls (LFC = 0.788±0.3317, 
— 1.73 fold, jft = 0.035) (Fig. 3), perhaps not surprising considering 
this value in the microarray experiment for the RA was ~ 1 fold 
(LFC — 0.555). No significant changes in expression of Thrb (Fig. 4) 
Gygl (Fig. 5) and Pdk4 (Fig. 6) were observed between the BPA 
exposed and control groups. Collectively, these qRT-PCR results 
largely confirm the microarray expression changes. 




Figure 2. Altered fast myosin heavy chain expression. Real-time 
PCR quantification of myosin heavy chain 6, cardiac isoform alpha 
(IVlyh6) expression in the left ventricle (LV), right ventricle (RV), left 
atrium (LA) and the right atrium (RA) of rhesus monkey (Macaca 
mulatta) fetuses at birth. Fetuses were exposed maternally to a 400 |ig/ 
Kg. body vj/eight of BPA dose (B) or 150 ^L ethanol (C) during late 
gestation (LG, days 100±2-term). Data were analyzed by tissue and 
bars represent the mean ± S.E.IVl. of the log2 (fold change). A two 
sample t-test was performed to identify significant BPA effects and a 
two-tailed psO.05 was considered statistically significant. Marginal 
significance between the groups was indicated with a '#' symbol. 
doi:10.1371/journal.pone.0089096.g002 

representing approximately 7% of the total myosin protein in 
human adult left ventricles [38-40]. Importantly, down-regulation 
of Myh6 to nearly undetectable levels in the failing adult human 
left ventricles has been noted [38,40], indicating that this myosin 
isoform influences ventricular myocardial fitness. Research on 
breast cancer cell lines indicates an inhibitory effect of 17fi- 
estradiol, BPA and other known endocrine disruptor compounds; 
genistein and polychlorinated biphenyl congener 54 (PCB54), on 
the Mjh6 transcript levels [41]. Although, the functional signifi- 
cance of Myh6 in breast cancer cells is unknown, these studies do 



Tissue glycogen content 

Glycogen is an important source of glucose and hence energy in 
the newborn fetus, therefore we measured the tissue glycogen 
content towards gaining insight into the mechanism of Mjh6 
down-regulation at the metabolite level. However, ventricular 
muscle glycogen, in the LG maternal BPA exposed vs. matched 
control fetuses was not significandy different (Fig. 7). 

Discussion 

We surveyed the broad-spectrum effects of maternal BPA 
exposure on female fetal cardiac gene expression profiles in a non- 
human primate using microarray and qRT-PCR analyses. We 
report evidence of BPA-induced gene expression changes in the 
primate fetal cardiovascular tissue indicating potential altered 
cardiac development. This study provides insights into the 
potential consequences of BPA exposure on human development 
and long-term health. 

Maternal BPA exposure during LG causes a down-regulation of 
fast myosin heavy chain isoform of mammalian myocardium 
{Mjh6) mRNA expression in the LV of rhesus fetal heart. Myh6 
gene products' contribution to total myocardial myosin, vary by 
tissue (ventricle vs. atrium), species (rodents vs. humans) and 
developmental stage (fetal vs. adult) [37-39]. While less abundant 
in the fetal ventricles, Mjh6 induction normally occurs just prior to 
birth and predominates in the rodent ventricles (>95%) while 
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Figure 3. Expression changes in Adam12-I mRNA. Real-time PCR 
quantification of the long isoform of A Disintegrin And Metalloprotease 
(Adam12-I) expression in the left ventricle (LV), right ventricle (RV), left 
atrium (LA) and the right atrium (RA) of rhesus monkey (Macaca 
mulatta) fetuses at birth. Fetuses were exposed maternally to a 400 |ig/ 
Kg. body weight of BPA dose (B) or 150 ^L ethanol (C) during late 
gestation (LG, days 100±2-term). Data were analyzed by tissue and 
bars represent the mean ± S.E.M. of the log2 (fold change). A two 
sample t-test was performed to identify significant BPA effects and a 
two-tailed p£0.05 was considered statistically significant. Asterisk (*) 
indicate significantly different groups. 
doi:10.1371/journal.pone.0089096.g003 
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Figure 4. Real-time PCR quantification of thyroid Kiormone 
receptor beta (Tlirb) expression. Transcript levels of Thrb in the left 
ventricle (LV), right ventricle (RV), left atrium (LA) and the right atrium 
(RA) of rhesus monkey (Macaca mulatta) fetuses at birth. Fetuses were 
exposed maternally to a 400 |ig/Kg. body weight of BPA dose (B) or 
150 |j,L ethanol (C) during late gestation (LG, days 100±2-term). Data 
were analyzed by tissue and bars represent the mean ± S.E.M. of the 
log2 (fold change). A two sample t-test was performed to identify 
significant BPA effects and a two-tailed psO.05 was considered 
statistically significant. 
doi:10.1371/journal.pone.0089096.g004 

provide a model of how BPA may down-regulate Myh6. Given that 
neonatal cardiomyocytes express functional estrogen receptors 
[23] it is plausible that hormone mediated pathways drive the BPA 
induced inhibition oi Myh6 expression observed in our study. 

Thyroid hormone status is an important determinant of myosin 
isoform composition in the cardiac ventricles [42]. Initiation of 
thyroid hormone receptor alpha binding of Tg at birth initiates the 
onset of postnatal cardiac phenotype [43] and myosin heavy chain 
isoforms in postnatal ventricles are direct transcriptional targets of 
thyroid hormone as indicated by reduced Myh6 expression in 
hypothyroid adult ventricular myocardium [44-47]. Research in 
rodent cells revealed that BPA acts as a thyroid hormone 
antagonist and represses thyroid hormone transcriptional target 
genes by non-competitively binding to the thyroid hormone 
receptor, and recruiting a co-repressor, N-CoR [6,7]. Our result 
supports this mechanism of action and indicates that BPA inhibits 
the induction oi Myh6 and may interfere with the transformation 
of the fetal heart into postnatal heart. Another non-exclusive 
mechanism of Myh6 down-regulation may be the repression of 
thyroid hormone receptor expression [48]. However, we did not 
observe any significant differences in the mRNA expression of 
these receptors in response to BPA exposure. A recent report 
showed that miR-205 silences TRAP220/MED1, a co-activator of 
several nuclear hormone receptors including thyroid hormone 
receptor, in hypoxic human placental trophoblasts [49] . Induction 
of miR-205 in our study is particularly important as TRAP220 is 
indispensable for embryonic development and homozygous 
TRAP220 knockout mice experience placental insufficiency and 
defective heart development, consistent with heart failure, and 
suffer early embryonic mortality [50,51]. In vitro studies conducted 
on embryonic fibroblasts from these TRAP220 knockout mice 
indicated perturbed thyroid hormone signaling. The BPA 
mediated induction oi miR-205 and miR-224 in fetal cardiac tissue 
is compelling as these are known to regulate tissue remodeling 
[49,52]. 

Maternal BPA exposure during LG up-regulates Adaml2-l 
expression in the LV (microarray), RV and RA (microarray and 



Figure 5. Real-time PCR quantification of glycogenin (Gyg1) 
expression. Expression of Gygl mRNA in the left ventricle (LV), right 
ventricle (RV), left atrium (LA) and the right atrium (RA) of rhesus 
monkey (Macaca mulatta) fetuses at birth was measured. Fetuses were 
exposed maternally to a 400 |jg/Kg. body weight of BPA dose (B) or 
150 jiL ethanol (C) during late gestation (LG, days 100±2-term). Data 
were analyzed by tissue and bars represent the mean ± S.E.M. of the 
log2 (fold change). A two sample t-test was performed to identify 
significant BPA effects and a two-tailed psO.05 was considered 
statistically significant. 
doi:10.1371/journal.pone.0089096.g005 

qRT-PCR) of rhesus monkey fetuses. Adaml2 is a member of the 
ADAMs multi-gene family that is expressed during early 
embryonic and fetal stages of development and in the adult heart 
[53,54] with demonstrated roles in the integration of ceU-to-ceU 
and ceU-to-extraceUular matrix interaction, as well as growth 
factor signaling [55]. Adaml2 has been implicated in the 
development of cardiac hypertrophy [56] and elevated expression 
of Adaml2 is observed in human and rodent cardiac tissue during 
hypertrophic obstructive cardiomyopathy [53] and agonist- 
induced cardiac hypertrophy [57,58]. Although the physiological 
pathway associated with up-regulation oiAdaml2-l is not clear, we 
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Figure 6. Pyruvate dehydrogenase kinase isozyme 4 (Pdk4) 
mRNA expression. Real-time PCR quantification of Pdk4 gene was 
performed on the left ventricle (LV), right ventricle (RV), left atrium (LA) 
and the right atrium (RA) samples of rhesus monkey (Macaca mulatta) 
fetuses at birth. Fetuses were exposed maternally to a 400 |ig/Kg. body 
weight of BPA dose (B) or 150 |iL ethanol (C) during late gestation (LG, 
days 100±2-term). Data were analyzed by tissue and bars represent the 
mean ± S.E.M. of the log2 (fold change). A two sample t-test was 
performed to identify significant BPA effects and a two-tailed psO.05 
was considered statistically significant. 
doi:10.1371/journal.pone.0089096.g006 
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Figure 7. Tissue glycogen estimation. Ventricular glycogen content 
was assessed using a calorimetric method on the left ventricle (LV) and 
the right ventricle (RV) of rhesus monkey (Macaca mulatta) fetuses at 
birth. Fetuses were exposed maternally to a 400 |ig/Kg. body weight of 
BPA dose (B) or 150 fiL ethanol (C) during late gestation (LG, days 100- 
term). Data were analyzed by tissue and bars represent the mean ± 
S.E.M., of the tissue glycogen content measured in ^g for every 
milligram tissue. A two sample t-test at a two-tailed psO.05 was 
considered statistically significant. 
doi:1 0.1 371 /journal.pone.0089096.g007 

speculate that the transforming growth factor beta (TGF(3) 
pathway may be involved. Research conducted using hepatic 
stellate cells [59,60], fibroblasts and mammary epithelial cells [61], 
reported that TGF[3, isoform 1 (TGF(3l) is a transcriptional 
activator of both the long and the short isoforms o{ Adaml2 that 
reciprocally promotes TGFp/smad signaling by enhancing the 
TGF(3 type II receptor internalization [62]. Further, BPA up- 
regulates the expression of other TGFfi isoforms, TGFp2 and 
TGFfiS transcripts but not TGFpl, in neural tissue [63] and 
endometrial cancer cells [64], respectively. Additionally, up- 
regulation of miR-224 observed in our study further supports this 
interpretation as TGF(31 was shown to activate miR-224 
transcription in mouse granulosa cells [52]. Lack of any 
significantly increased TGFfi mRNA expression in our microarray 
analysis (data not shown) is however not surprising as temporal 
differences may exist in the expression of both the Adaml2-l and 
TGFj] transcripts. 

Our microarray analysis shows that BPA caused largely distinct 
gene expression changes in different cardiac quadrants and at both 
the developmental stages. Indeed, different quadrants of the 
mammalian heart (LV, RV, LA & RA) are divergent in structure, 
contractile function and metabolism as confirmed by compart- 
ment specific myocardial gene expression profdes [65-67]. 
Moreover, cardiac remodeling events that underlie important 
cardiovascular pathophysiologies vary in accordance with the 
specific heart compartments [67-69]. Hence, the influence of BPA 
on different anatomical subdivisions of the heart is presumed to be 
variable. Additionally, we observed larger 'within group' variabil- 
ity in gene expression, which may signify the genetic heterogeneity 
in the sample groups studied. While gene expression studies would 
benefit from a larger sample size, nonetheless our findings offer the 
first evidence of the potential of the maternal BPA exposure to 
alter cardiac development and health in non-human primate 
rhesus monkey. 

In conclusion, our data indicate that chronic exposure to 
maternally delivered BPA significandy alters the transcript 
expression profiles of several coding and non-coding genes in 
rhesus monkey fetal heart tissue. The results of this study are 
important for multiple reasons. First, the observed changes of 
Myh6 and Adaml2-l mRNA expression are consistent with events 



that underlie cardiac remodeling during cardiovascular patho- 
physiologies such as cardiac hypertrophy [38,40,56]. Secondly, 
similarities in the metabolism of rhesus monkey and humans make 
these observations relevant to the latter, and may provide 
physiological insights into the recendy reported correlations 
between urinary BPA concentration and the incidence of 
cardiovascular disease [20,21]. Finally, our results emphasize the 
impact of physiologically relevant BPA exposure as the average 
maternal serum unconjugated BPA level in our study was 
0.68 ng/ml [29], which is closer to the levels reported in human 
maternal serum during gestation [70,71]. The results of this study 
conflict with recent reports that suggest that BPA poses no human 
health risk in either adults or newborns [72]. This research adds to 
the evidence suggesting that the rising incidence of cardiovascular 
and metabolic diseases in the human population may be in part a 
reflection of our lifetime dietary and environmental exposure to 
chemicals such as BPA. 
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